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Abstract 

Motivated by the possibility of comparing theoretical predictions of Lake Vostok's composition with future in situ 
measurements, we investigate the composition of clathrates that are expected to form in this environment from the air 
supplied to the lake by melting ice. In order to establish the best possible correlation between the lake water com- 
position with that of air clathrates formed in situ, we use a statistical thermodynamic model based on the description 
of the guest-clathrate interaction by a spherically averaged Kihara potential with a nominal set of potential parame- 
ters. We determine the fugacities of the different volatiles present in the lake by defining a "pseudo" pure substance 
dissolved in water owning the average properties of the mixture and by using the Redlich-Kwong equation of state 
to mimic its thermodynamic behavior. Irrespective of the clathrate structure considered in our model, we find that 
xenon and krypton are strongly impoverished in the lake water (a ratio in the 0.04-0. 1 range for xenon and a ratio 
in the ~0.15-0.3 range for krypton), compared to their atmospheric abundances. Argon and methane are also found 
depleted in the Lake Vostok water by factors in the 0.5-0.95 and 0.3-0.5 ranges respectively, compared to their at- 
mospheric abundances. On the other hand, the carbone dioxide abundance is found substantially enriched in the lake 
water compared to its atmospheric abundance (by a factor in the 1.6-5 range at 200 residence times). The comparison 
of our predictions of the CO2 and CH4 mole fractions in Lake Vostok with future in situ measurements will allow 
disentangling between the possible supply sources. 
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1. Introduction 

Over 100 liquid water lakes have been identified be- 
neath the ice sheets of Antarctica (Smith et al. 2009). 
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4 Among them, Lake Vostok lying buried beneath ~4 

5 kilometers of ice is the largest known subglacial lake 
e on Earth. Estimates for the age of this lake range from 
7 one million years (Kapitsa et al. 1996) to 15 million 
a years (Bell et al. 2002). Because Lake Vostok is under 
9 a very thick ice sheet, it is unlikely that gases within the 
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lake water are in equilibrium with air. From a biologi- 
cal point of view, the study of Lake Vostok is interesting 
because it provides an environment that has been sealed 
off from light and atmosphere for possibly several mil- 
lions of years. As such, the Lake Vostok environment 
could be analogous to that of the Jovian icy moon Eu- 
ropa, which is envisaged to harbor a subsurface ocean 
with conditions compatible with habitability (Hand et 
al. 2006; Pasek & Greenberg 2012). 

After cumulated decades of effort, researchers have 
succeeded in drilling through 4 kilometers of ice to the 
surface of the subglacial Lake Vostok on 5 February 
2012 (Showstack 2012). The precise sampling of Lake 
Vostok in a near future will bring important constraints 
on the delivery processes of volatiles to the lake and the 
chemical pathways that may, or may not make life pos- 
sible deep below the ice sheet. 

Ice overlying the lake is characterized by air bubbles 
trapped at the air/ice interface that are compressed with 
increasing depth. As a consequence, these air bubbles 
transform to clathrate hydrates (hereafter clathrates) be- 
low a certain depth, depending on the temperature of the 
ice sheet. In the contact region between the bottom of 
the ice sheet and the lake surface, most of the air in the 
ice sheet is trapped within these clathrates. In situ mea- 
surements have shown that the average N2/O2 ratio of 
the clathrates in this contact region approaches the at- 
mospheric value (Ikeda et al. 1999). Melting of these 
clathrates is thus responsible for a transfer of most of 
the atmospheric air to the Lake Vostok water through 
the ice sheet (Lipenkov & Istomin 2001). However, in 
the contact region, clathrates could also be stable in the 
lake water thus reducing the amount of gases in aque- 
ous solution. This might have an impact on the oxygen 
availability for possible life in the lake and it is thus very 
important to thoroughly investigate the composition of 
Lake Vostok. 

The lake composition was previously investigated by 
McKay et al. (2003) (hereafter MK03). Based i) on 
the prediction that air clathrate is stable at depth higher 
than ~1500 m in sub-glacial lakes (Lipenkov & Istomin 
2001) and ii) on the fact that selective gas trapping oc- 
curs during the formation of clathrates (van der Waals 
& Platteeuw 1959; Miller 1974; Lunine & Stevenson 
1985), MK03 estimated the level of gases dissolved in 
Lake Vostok as a function of its evolution. To do so, they 
modeled the balance between the gas (i.e. air) supplied 
by ice melting and gas trapped by clathrate, and assum- 
ing that the lake is a closed system, with no gas escape. 
However, the approach followed by MK03 to calculate 
the relative partitioning of volatiles trapped in clathrate 
is not adapted for a gas mixture in which more than two 



species are included (Lunine & Stevenson 1985), as it 
is the case for air clathrate. Also, the formalism used 
by the authors is a crude approximation of the statistical 
thermodynamic model that is usually used to investigate 
the composition of multiple guest clathrates. Indeed, 
it is based on the assumption that the Langmuir con- 



stants (see Section 2.2 1 are very close to each other for 
the different gases, irrespective of the considered type 
of clathrate cages (Lunine & Stevenson 1985). 

In this paper, motivated by the possibility of compar- 
ing theoretical predictions of Lake Vostok composition 
with in situ measurements in the near future, we reinves- 
tigate the assumptions of MK03 by considering a more 
complete set of species dissolved in the lake. We thus 
determine more accurately the composition of clathrates 
that should form in this environment. In order to estab- 
lish a better correlation between the lake composition 
with that of air clathrates formed in situ, we use a statis- 
tical thermodynamic model based on experimental data 
and derived from the original work of van der Waals 
& Platteeuw (1959). This approach, used today in in- 
dustry and science, has saved substantial experimental 
effort for the determination of: (i) the equilibrium pres- 
sure of a clathrate formed from various mixtures; and 
(ii) the mole fraction of the different species trapped in 
the clathrate from a given fluid phase (see Sloan & Koh 
(2008) for details). The major ingredient of our model 
is the description of the guest-clathrate interaction by 
a spherically averaged Kihara potential with a nomi- 
nal set of potential parameters. Section [2] is devoted to 
the description of i) the strategy utilized to compute the 
lake composition and ii) the statistical thermodynamic 
model used in this work. In Section [3] we present the 
results of our calculations concerning the time evolu- 
tion of clathrate and lake compositions. In Section [4j 
we compare our results with those of MK03 and dis- 
cuss their implications for the physical and biological 
processes that may take place in Lake Vostok and anal- 
ogous environments such as the putative internal ocean 
of the Jovian icy moon Europa. 

2. Modeling the lake composition 

2.1. Computational approach 

Our basic assumptions are similar to those formulated 
by MK03. We assume that Lake Vostok is a closed sys- 
tem at constant pressure of 35 MPa with a mean tem- 
perature of -3°C, and is physiographically stable over 
time (Lipenkov & Istomin 2001) for a description of 
Lake Vostok's thermodynamic conditions). In our sys- 
tem, water and air are delivered to the lake when melting 



occurs due to the slow downward motion of the over- 
lying gas-rich ice layers, and gas-free water leaves the 
lake as ice accretes to the bottom of the ice sheet in 
regions where ice moves outward (Jouzel et al. 1999; 
Bell et al. 2002; MK03). Air is assumed to be supplied 
to the lake at a concentration of about 90 cm 3 at STP 
(T = 273.15 K and P = 0.1013 MPa) per kg of melted 
ice (Jouzel et al. 1999; Lipenkov & Istomin 2001). We 
express the age of the lake in units of the residence time 
(denoted by RT in the following), defined as the time 
required for all water of the lake to be renewed through 
the melting-freezing cycle of water flowing through the 
volume of the lake. Estimates of the RT for Lake Vostok 
range from 5,000 years (Philippe et al. 2001) to 125,000 
years (Siegert et al. 2003). 

Because the clathrate formed in Lake Vostok is crys- 
tallized from dissolved air, it is expected to be domi- 
nated by N2. A good approximation of its dissociation 
pressure is to assume it is equivalent to that of pure N2 
clathrate. The dissociation pressure P~l ss of N2 clathrate 
follows an Arrhenius law (Miller 1961) defined as 



log{Pt S ) = A + 



B 



(1) 



132 where P^ s is expressed in Pa and T is the temperature 

133 in K. The constants A and B are set to 9.86 and -728.58, 

134 respectively, and fit to experimental data (Thomas et al. 

135 2007). At 270. 15 K, P d £ s is found to be -14.6 MPa and 

136 the corresponding fugacity is ~ 13.7 MPa (see Sec. |2.3[ ). 

137 However, at high hydrostatic pressure (here 35 MPa), 

138 the calculated dissociation fugacity needs a correction 

139 due to the specific volume difference between clathrate 

140 and the liquid solution. The effective dissociation fu- 

141 gacity of clathrate is found equal to ~21.2 MPa (see Eq. 

142 [T3"| l, which translates into an equivalent gas pressure of 

143 ~22.5 MPa. The pressure of dissolved N2 gas must in- 

144 crease with time up to this value to enable the trigger- 

145 ing of clathrate formation in Lake Vostok. After having 

146 computed the Henry's constant of N2 at the hydrostatic 

147 pressure of 35 MPa (see Sec. |2.3| >, we find that the gas 
us pressure of 22.5 MPa corresponds to a volume of ~ 1735 

149 cm 3 (STP) of N2 per kg of water. These estimates im- 

150 ply that, at ~25 RT, the concentration of air dissolved 

151 in Lake Vostok reaches a solubility limit of ~2220 cm 3 

152 (STP) per kg of water, a value lower than the one found 

153 by MK03 (2500 cm 3 at STP), just before the beginning 

154 of clathrate formation. 

155 Since the time for dissolved gases to diffuse through 

156 the lake is smaller than 1 RT, and due to the long 

157 timescales involved (several tens of RT, see hereafter), 

158 the gas may be supposed to be uniformly mixed in the 

159 lake water at 25 RT. Thereafter, the total amount of dis- 
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Figure 1 : Structural scheme of our computational approach. 

solved gas in Lake Vostok remains approximately con- 
stant and excess gas continuously forms clathrate. Be- 
cause the composition of clathrate may depart from that 
of the coexisting gas phase, progressive clathrate forma- 
tion could in turn influence the composition of the gas 
dissolved in Lake Vostok. Table[T]gives the composition 
of the air dissolved in the lake water at the beginning of 
our computation (i.e. at RT). 

Based on this approach, we have elaborated a compu- 
tational procedure (see Figure 1) aiming at calculating 
the volume of each gas dissolved in Lake Vostok (per 
unit of kg of lake water) and trapped in the newly form- 
ing air clathrates, as a function of the temporal evolution 
of the lake (in units of RT). The procedure is depicted 
as follows: 



L. 



At time t — 0, the mole fraction x£ (0) of species 
K, relative to the sum of all mole fractions of gases 
dissolved in the lake water, is given in Table[T] The 
corresponding volume of gas V l £ ke {0) per unit of kg 
of water is determined. The mole fraction x c ^ ath (Q) 
and the volume V^ a>h (Q) of the enclathrated species 
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Table 1: Gas dissolved in Lake Vostok. From left to right: composition of dry air by volume (Lodders & Fegley 1998), critical temperatures and 
pressures (Lide 2002), partial molar volumes ("Rudakov et al. 1996; * Anderson 2002; ' Kennan & Pollack 1990), Henry's constants at STP (Sander 
1999) and calculated at 35 MPa of hydrostatic pressure and 270.15 K, and fugacities calculated at 25 RT and at pressures of 1 atm and 35 MPa, 
respectively. Properties of species that are not incorporated in clathrate are intentionally left blank (see text). 



Species K 


Y air 
X K 


T c 




Vk 


Hlatm,K 


Hkydro,K 


f\ulm,K 




fhidro,K 








(K) 


(MPa) 


(cm 3 /mol) 


(L atm mol" 1 ) 


(L atm mol -1 ) 


(MPa) 




(MPa) 




N 2 


0.7808 


126.19 


3.40 


36" 


1639.34 


2868.19 


11.80 




18.29 




o 2 


0.2095 


154.48 


5.04 


33" 


769.23 


1284.54 


1.49 




2.30 




Ar 


9.34 x 10~ 3 


150.87 


4.90 


33" 


714.29 


1192.79 


6.15 x 10 


-2 


9.53 x 10- 


-2 


co 2 


3.50 x 10" 4 


304.25 


7.38 


38.4* 


28.57 


51.88 


9.22 x 10 


-5 


1.43 x 10- 


-4 


Ne 


18.18 x 10~ 6 




















He 


5.24 x 10~ 6 




















CH 4 


1.70 x 10- 6 


190.6 


4.6 


35" 


714.29 


1230.44 


1.12 x 10 


-5 


1.73 x 10- 


-5 


Kr 


1.14 x 10- 6 


209.4 


5.5 


30" 


400 


637.52 


4.20 x 10 


-6 


6.51 x 10- 


-6 


H 2 


0.55 x 10- 6 




















CO 


1.25 x 10~ 7 


133.1 


3.5 


37" 


1010.10 


1794.95 


1.16 x 10 


-6 


1.80 x 10- 


-6 


Xe 


0.87 x 10~ 7 


289.77 


5.84 


AT 


232.56 


482.72 


1.86 x 10 


-7 


2.89 x 10- 


-7 
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K are set to zero; 

at t = t + 1, a subroutine (depicted in Section 2.2 1 
allows us to compute the mole fraction x c Jf th (t) of 
each enclathrated guest K by using as input its vol- 
ume (per unit of kg of water) in Lake Vostok at 
time t - 1 ; 

the new volumes V C J? if) and V^ ke (t) are calcu- 
lated for each species K. The new mole fraction 
x l £ ke (t) is obtained and normalized before being 
reintroduced in the next loop, until time t reaches 
150 RT, a value considered as infinity time. In the 
following, we consider that infinity time is the time 
at which the lake composition reaches steady state. 



2.2. The statistical thermodynamic model 

223 

To calculate the relative abundances of guest species 224 

incorporated in a clathrate from a coexisting gas of SpeC- 225 

ified composition at given temperature and pressure, our 220 

subroutine follows the method described by Lunine & 227 

Stevenson (1985), Thomas et al. (2007, 2008, 2009) 228 

and Mousis et al. (2010), which uses classical statisti- 229 

cal mechanics to relate the macroscopic thermodynamic 230 
properties of clathrates to the molecular structure and 
interaction energies. It is based on the original ideas 
of van der Waals & Platteeuw (1959) for clathrate for- 
mation, which assume that trapping of guest molecules 

into cages corresponds to the three-dimensional gener- 231 

alization of ideal localized adsorption. This approach is 232 

based on four key assumptions: 233 

234 

1 . The host molecules contribution to the free energy 235 

is independent of the clathrate occupancy. This as- 230 

sumption implies that the guest species do not dis- 237 

tort the cages; 23s 



2. the cages are singly occupied and guest molecules 
rotate freely within the cage; 

3. guest molecules do not interact with each other; 

4. classical statistics is valid, i.e., quantum effects are 
negligible. 

In this formalism, the fractional occupancy of a guest 
molecule K for a given type q (q = small or large) of 
cage in a clathrate structure can be written as 



yK,q = 



CK,q fhydro,K 
1 + 2/ Cj^ q f h y, 



(2) 



dro,J 



where the sum in the denominator includes all the 
species that are present in the gas dissolved in lake wa- 
ter. Cx, q is the Langmuir constant of species K in the 
cage of type q, and fh y dro,K me fugacity of guest species 
K computed at the considered hydrostatic pressure (see 
Sec. [231). 

The Langmuir constant Cx,q depends on the strength 
of the interaction between each guest species and each 
type of cage, and can be determined by integrating the 
molecular potential energy within the cavity as 



An 



(3) 



where R c represents the radius of the cavity assumed 
to be spherical, kg the Boltzmann constant and w^ 9 (r) 
is the spherically averaged potential (here Kihara or 
Lennard-Jones potential) representing the interactions 
between the guest molecules K and the H2O molecules 
forming the surrounding cage q. This potential w(r) can 
be written for a spherical guest molecule, as (McKoy & 
Sinanoglu 1963) 
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Table 2: Parameters for the clathrate cavities. R c is the radius of the 
cavity (values taken from Parrish & Prausnitz 1972). b represents 
the number of small (b s ) or large (b[) cages per unit cell for a given 
structure of clathrate (I or II) with volume V c , z is the coordination 
number in a cavity. 



Clathrate structure 


I 


II 


N w 


46 


136 


V c (A 3 ) 


12.0 3 


17.3 3 


Cavity type 


small large 


small large 


Rc (A) 


3.975 4.300 


3.910 4.730 


b 


2 6 


16 8 


z 


20 24 


20 28 



269 
270 



12 



(4) 274 

239 with the mathematical function 5 N (r) in the form 



'OK-sr-K-sn-™ 

In Eqs. (|4]) and PJ, z is the coordination number of the 
cell. 

Parameters z and R c depend on the type of the cage 
(small or large) and also on the structure of the clathrate. 
Indeed, the air clathrates considered here may have 
structure I or II that are characterized by different ar- 
rangement of small and large cages of different sizes 
(Sloan & Koh, 2008). The corresponding parameters 
z and R c are given in Table [2] The intermolecular pa- 
rameters a, <j and e describing the guest molecule-water 
(K— W) interactions in the form of a Kihara or Lennard- 
Jones potential are listed in Table [3] 

Finally, the mole fraction x K ar of a guest molecule K 
in a clathrate of a given structure (I or II) can be calcu- 
lated with respect to the whole set of species considered 
in the system as 



276 
277 
278 
279 



b s y K , s + b e y K ,e 

bs'Ljyis + beTijyj/ 



(6) 



where b s and b/ are the number of small and large cages 
per unit cell, respectively, for the clathrate structure un- 
der consideration. Note that the sum of the mole frac- 
tions of the enclathrated species is normalized to 1 . We 
also neglect Ne, He and H2 in our computation of the 
composition of clathrates as these species have low trap- 
ping propensities in these structures (Lunine & Steven- 
son 1985; Sloan & Koh 2008). 



2.3. Determination of the fugacities 

To determine fhydro,K< we first calculate the specific 
volume v of the considered mixture via the Redlich- 
Kwong equation of state (Redlich & Kwong 1949): 



P = 



R T 



b y/fv(v + b) 



(7) 



with 



R 2 r 2 5 

0.42748 — and b = 0.08664 



R T r 



(8) 



where R is the gas constant, T the ambient tempera- 
ture, T c and P c the critical temperature and pressure of 
the substance (see Table [TJ, and P its vapor pressure. 
Just as any other cubic equation of state, the Redlich- 
Kwong EOS has to be applied only to pure substances. 
For mixtures, however, the same equation is applied but 
certain mixing rules are applied to obtain parameters a 
and b, which are functions of the properties of the pure 
components. This corresponds to the creation of a new 
"pseudo" pure substance that has the average properties 
of the mixture. Following Redlich & Kwong (1949), we 
use the mixing rules: 



a m = EiZ; x' flfe x l f e ay with a u 



bin — %i b{ 



(9) 



Here the mixture pressure P„, dissolved in water cor- 
responds to the sum of the individual gas pressures Pk- 
Each Pk is expressed as a function of the Henry's law 
coefficient as: 



tt „ lake 

tlhydro,K X y K , 



(10) 



where Hhydw,K is the Henry's law constant calculated 
for species K at the lake's hydrostatic pressure and am- 
bient temperature, and y l ^ ke the mole fraction of species 
K expressed per unit of volume of water. This term is 
derived from (Krichevsky & Kasarnovsky 1935): 



In 



H, 



hydro, K 



H 



iatm.K 



^(Phydro- 1.013 X10 5 ), 



(ID 



where Phydro is the lake's mean hydrostatic pressure 
(~35 MPa), Vk the partial molar volume of species K 
and H, 

am,K the Henry's constant determined at 1 atmo- 
sphere and at 0°C (see Table [TJ. 

The fugacity of the mixture at 1 atm pressure f\ a tm,m 
is then determined from the following relation (Redlich 
& Kwong 1949): 



5 



Table 3: Parameters for Kihara and Lennard-Jones potentials. <tk-w is the Lennard-Jones diameter, £k-w is the depth of the potential well, and 
aic-w is the radius of the impenetrable core, for the guest-water pairs. 



Molecule 


a- K -w (A) 


e K-w/k B (K) 


a K - W (A) 


Reference 


co 2 


2.97638 


175.405 


0.6805 


Sloan & Koh (2008) 


N 2 


3.13512 


127.426 


0.3526 


Sloan & Koh (2008) 


o 2 


2.7673 


166.37 


0.3600 


Parrish & Prausnitz (1972) 


CH 4 


3.14393 


155.593 


0.3834 


Sloan & Koh (2008) 


CO 


3.1515 


133.61 


0.3976 


Mohammadi et al. (2005) 


Ar 


2.9434 


170.50 


0.184 


Parrish & Prausnitz (1972) 


Kr 


2.9739 


198.34 


0.230 


Parrish & Prausnitz (1972) 


Xe 


3.32968 


193.708 


0.2357 


Sloan & Koh (2008) 
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flatm.m ~ 

: 

dm I 1 



RTV 2 \v m + h 



+ In 
In 



RT 



1 , (Vm + b m ) 



(12) 



297 Because it is calculated at the hydrostatic pressure, 

298 fhydro,m 

is related to f\ atm , m via (Miller 1974): 



In 



fhydro,m 



p hydro -\m r i x io : 

R T 



flatm.m 
,5 



= — n In aw 



■(" Vh 2 ~ Vc,m), 



(13) 



where aw is the activity of water relative to pure liquid 
water (here aw is ~1), n is the moles of water per mole 
of species K in the clathrate (n ~6), Vh 2 o ls me molar 
volume of liquid water (1.8 x 10~ 5 m 3 ) and V cjn is the 
volume of clathrate that contains 1 mole of substance m 
(here approximated to be dominated by N2). Finally the 
fugacity coefficient cf> of the mixture, which is defined 
as the ratio of fh y dro,m to Phydro,m> allows to retrieve the 
individual fugacities following 
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fhydro,K 



K- 



(14) 



3. Results 



3.1. Evolution of clathrate composition 3 „ 

The composition of structure I and II clathrates that 34s 

might form from the air dissolved in Lake Vostok has 349 

been computed using the approach presented above. If 350 

the lake has experienced less than ~25 RT, the differ- 351 

ent species remain completely in solution. At longer 352 

times, clathrates start to form. Figures 2 and 3 display 353 

the compositions of the clathrates of same structure sue- 354 

cessively formed at different epochs of the lake evolu- 355 

tion. At each timestep, the total volume of enclathrated 356 



gases is 90 cm 3 . Irrespective of the clathrate struc- 
ture, the mole fractions of volatiles trapped in clathrates 
converge toward their atmospheric mole fractions at in- 
finite RT as the result of mass conservation (MK03). 
Steady state is reached 40-80 RT after clathrate for- 
mation. In both structures, N 2 and 2 remain the 
main gases trapped in clathrates with variations of their 
mole fractions depending on the considered structure 
and epoch of clathrate formation. In the case of struc- 
ture I clathrates, the mole fractions of trapped N 2 and 
2 change slightly over time. Indeed, the mixing ra- 
tio of N 2 decreases from ~0.86 in clathrates formed 
at early epochs to 0.78 in clathrates formed at infinite 
time whereas that of 2 simultaneously increases from 
~0.12 to 0.21 in clathrates successively formed dur- 
ing the same time interval. In the case of Structure II 
clathrates, the mole fractions of N 2 and 2 incorporated 
in clathrates are even less sensitive to their formation 
epoch in the lake. Thus, the mole fraction of N 2 in- 
creases from 0.74 to 0.78 while that of 2 decreases 
from 0.24 to 0.21 in clathrates formed between first for- 
mation and steady state. 

Moreover, we note that, with an enrichment factor 
of ~22 compared to the atmospheric abundance (see 
Table [TJ, Xe is efficiently sequestrated at early forma- 
tion epochs (~25 RT) in structure I clathrates. To a 
lesser extent, a similar trend is calculated for Kr and 
CH4 whose mole fractions are both enriched by a factor 
of ~3.2 at early epochs of structure I clathrate forma- 
tion in the lake. The same remarks apply for the trap- 
ping of volatiles in structure II clathrates, for which the 
mole fractions of Xe, Kr and CH4 are enriched by fac- 
tors ~78, 6.7 and 2 at early formation epochs. On the 
other hand, the trapping efficiencies of Ar and CO in 
clathrate depend on the considered structure. In struc- 
ture I, the mole fraction of Ar remains similar to its at- 
mospheric value, irrespective of the considered RT. In 
structure II, however, we note that Ar is initially en- 
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Figure 2: Mole fractions of N 2 , 2 , Ar, C0 2 , CH 4 , Kr, Xe and CO 
trapped in structure I clathrates formed at different ages of the lake, 
with age expressed in units of RT. Clathrate formation starts at ~25 
RT (see text). 
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Figure 3: Mole fractions of N 2 , 2 , Ar, C0 2 , CH4, Kr, Xe and CO 
trapped in structure II clathrates formed at different ages of the lake, 
with age expressed in units of RT. Clathrate formation starts at ~25 
RT (see text). 
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riched by a factor of ~2 in clathrate. In addition, CO 
is initially slightly enriched by a factor of ~ 1.2 in struc- 
ture I clathrate while it is impoverished by a factor of 
~0.8 in structure II clathrate. In both structures and at 
early epochs of clathrate formation, the mole fraction of 
trapped CO2 is lower than its atmospheric value (factors 
of ~0.6 and 0.15 in structures I and II, respectively). In 
both cases, the mole fractions of all volatiles trapped in 
clathrate progressively converge with time towards their 
atmospheric values. 



3.2. Evolution of the lake composition 

Figures 4-7 illustrate the evolution of the composi- 
tion of the gas phase dissolved in Lake Vostok as a 
function of its RT in the cases of structures I and II 
clathrates. In both cases, the mole fractions of the dif- 
ferent volatiles relative to the sum of all mole fractions 
of dissolved gases and present in the lake linearly in- 
crease as long as its lifetime has not exceeded ~25 RT. 
After this epoch, these volatiles have reached their sol- 
ubility limit in the lake and clathrate forms from the gas 
supplied in excess by ice melting. The resulting volatile 
content of the lake at a given time is then controlled 
by the balance between the sequestration of gas by the 
forming clathrate and the air supplied by the melting 
ice. 

In the two cases, N2 and O2 remain the main gases 
dissolved in water. When considering the formation of 
structure I clathrate, and irrespective of the considered 
epoch, the N2 mole fraction in the lake water decreases 
from its atmospheric value (0.78) down to 0.66 at steady 
state. On the other hand, the O2 mole fraction slightly 
increases with time in water from 0.21 (the atmospheric 
mole fraction) to 0.33 in the lake. As a result, the N2/O2 
ratio in the lake water is much lower than the N2/O2 
ratio in structure I clathrate (see Fig. 6). Because the 
mole fractions of Xe, Kr, CH4 and CO are all enriched 
in structure I clathrate, they become correspondingly 
depleted by factors -0.04, 0.3, 0.3 and 0.8 in the lake 
water compared to their atmospheric abundances (see 
Table [TJ, and irrespective of the considered RT. In con- 
trast, the mole fraction of Ar in the lake water remains 
very close (by a factor of ~0.95) to its atmospheric value 
while that of CO2 rapidly increases with time up to 1 .6 
times its atmospheric mole fraction at 200 RT. 

In the case of structure II clathrate formation in Lake 
Vostok, the N2 mole fraction in the lake water increases 
from its atmospheric value (0.78) up to 0.82 at steady 
state. Inversely, the O2 mole fraction decreases in wa- 
ter from its atmospheric value (0.21) down to 0.18 at 
steady state. The resulting N2/O2 ratio is then signif- 
icantly higher in the lake water than the N2/O2 ratio 



1.0 
0.8 
I 0.6 

u 
ro 

I 0.4 



0.2 

0.0 

10" 2 
10" 3 
10" 4 
10" 5 
10' 6 
10' 7 
10' 8 
10" 9 fc 



Ar 



CO, 



■CO- 



■Xe- 



20 40 60 80 100 120 
Time (per unit of residence time) 

Figure 4: Mole fractions of N 2 , 2 , Ar, CO2, CH 4 , Kr, CO and Xe 
dissolved in Lake Vostok calculated in the case of structure I clathrate 
formation and expressed as a function of the age of lake, with age 
expressed in units of RT. 
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Figure 5: Mole fractions of N 2 , 2 , Ar, C0 2 , CH 4 , Kr, CO and Xe 
dissolved in Lake Vostok calculated in the case of structure II clathrate 
formation and expressed as a function of the age of lake, with age 
expressed in units of RT. 
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Figure 6: Top panel: evolution of the N2/O2 ratio in lake and in struc- 
ture I clathrate as a function of time. Bottom panel: evolution of the 
mole fractions of Ar, CCb, CH4, Kr, CO and Xe relative to their at- 
mospheric values in lake as a function of time in case of structure I 
clathrate formation. 
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in structure II clathrate, irrespective of the considered 
RT (see Fig. 7). The mole fractions of Xe, Kr, CH4, 
and Ar also become rapidly impoverished by the factors 
~0.1, 0.15, 0.5 and 0.5 in the lake water, respectively, 
compared to the atmospheric values (see Table [TJ, and 
irrespective of the considered RT. On the other hand, 
CO becomes moderately enriched by a factor ~1.3 in 
the water compared to its atmospheric value. Moreover, 
due to its lower propensity to be trapped in structure 
II clathrate, the CO2 mole fraction in the lake water 
presents the highest enrichment (up to a factor of ~5 
at 200 RT) compared to its atmospheric abundance. 



3.3. Clathrate density 



Once formed, clathrate may either float on the sur- 
face or sink to the lake floor according to its density 
with respect to that of the lake water. To investigate 
this point, we have computed the densities of structures 
I and II clathrates formed in Lake Vostok by following 
the method depicted in Sloan & Koh (2008; Eq. 5.2.1 
p268). We find that, in both cases, these clathrates are 
all lighter than liquid water (see Fig. 8). This conclu- 
sion is in agreement with the results obtained by MK03, 
although we do not calculate the same clathrate compo- 
sition. However, air clathrates have not been observed in 
the accreted ice above the lake (Siegert et al. 2000). To 
explain this, MK03 argued that significant amounts of 
CO2 could have been produced in Lake Vostok and sub- 
sequently incorporated in clathrates, thus changing their 
density with respect to that of clathrates simply formed 
in contact with air entering to the lake through the ice 
sheet. These clathrates would have sunk to the lake floor 
as a result of their higher density. We computed the min- 
imum fraction of CO2 present in the lake necessary to 
be trapped in air clathrates so that the clathrates sink to 
the lake floor. In the case of structure I clathrate, our 
model suggests that the minimum mole fraction of CO2 
present in Lake Vostok must be around ~0.05, with a 
corresponding mole fraction that is similar in clathrate 
(see Fig. 8). In the case of structure II clathrate, the min- 
imum mole fraction of CO2 needed to induce a higher 
clathrate density than that of lake water largely exceeds 
0.5, implying that this species becomes dominant in the 
lake. With such a high mole fraction of CO2 in the 
lake, structure I clathrate should be the most stable form 
(Sloan & Koh 2008) so we should expect a transition 
from structure II to structure I clathrate, which would 
present an even higher density. 
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Figure 8: Top panel: mole fraction of CO2 incorporated in structures I 
and II clathrates as a function of the mole fraction of CO2 dissolved in 490 
Lake Vostok. Bottom panel: densities of structures I and II clathrates 491 
as a function of the mole fraction of CO2 dissolved in Lake Vostok. 492 
The horizontal dashed line corresponds to the density of liquid water 
in Lake Vostok (Lipenkov & Istomin 2001). 
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The composition of Lake Vostok calculated with the 
present approach displays significant differences with 
the one predicted by MK03. These differences come 
from both the consideration of a larger set of molecules 
in the gas phase dissolved in the lake (we include CH4 
and CO) and from the use of a more sophisticated model 
to compute the mole fractions of all species simulta- 
neously trapped in the clathrates, in a self consistent 
way. As a consequence, our calculations suggest en- 
hancements of the Xe and Kr mole fractions trapped in 
clathrate by factors of ~ 7.6-26.9 and 2.Q-A.2, respec- 
tively, compared to the values previously determined by 
MK03. Our model therefore suggests that Xe and Kr 
abundances are ~0.04— 0.13 and 0.24-0.5 times those 
predicted by MK03 in Lake Vostok water. Our model 
also predicts an impoverishment of the Ar abundance 
in the lake water (factor ~0. 5-0.95 compared to the 
atmospheric abundance), in contrast with the enrich- 
ment (factor ~1.4) predicted by MK03. Moreover, our 
calculations show that CO2 is significantly enriched in 
the lake water (factor ~ 1.6-5 at 200 RT following the 
considered clathrate structure) compared to their atmo- 
spheric abundances. In both cases, CH4 appears de- 
pleted in the lake water compared to its atmospheric 
mole fraction (factor ~0.3-0.5). Depending on the con- 
sidered structure, CO can be found in the lake water 
moderately depleted (factor ~0.8 in the case of structure 
I clathrate) or enriched (factor ~ 1.3 in the case of struc- 
ture II clathrate) compared to its atmospheric value. It 
is important to note that, in absence of large amounts of 
CO2 dissolved in the liquid water, since the composition 
of the dissolved gas is dominated by N2 and because N2 
forms structure II clathrate (Sloan & Koh 2008), we ex- 
pect that structure II clathrate is the most likely to be 
formed in Lake Vostok. 

Note that the CO2 atmospheric abundance considered 
in this work corresponds to the one measured nowadays 
in the Earth's atmosphere, which is affected by indus- 
trial emissions. In the same way, the CO2 atmospheric 
abundance is estimated to have been 30% lower than the 
current value some 300 years ago and might have even 
been 10% higher 1-15 million years ago. Moreover, 
it has been inferred that extreme fractionation of gases 
could occur by formation of clathrates in Vostok ice due, 
at least partly, to different diffusion coefficients in ice of 
the considered gases (Ikeda et al. 1999). As a conse- 
quence, using the present atmospheric abundances for 
determining the composition of clathrates formed from 
air trapped for perhaps thousands of year in the ice sheet 
above the lake could be questionable. However, in situ 
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506 measurements have shown that the N2/O2 ratio deliv- 

507 ered to the lake at the bottom of the ice sheet is similar to 
5oa the present atmospheric ratio (Ikeda et al. 1999). Sim- 

509 ilarly, for CO2, it has been shown that the fractionation 

510 process in Vostok ice is smoothed out with depth, re- 

511 gaining the initial mean atmospheric concentration be- 

512 low the clathrate zone formation (Luthi et al. 2010). 

513 This information is unfortunately lacking for the other 

514 gases. However, additional tests performed by changing 

515 the CO2 composition in our approach show that, in any 

516 case, these variations do not affect our conclusions con- 

517 cerning the efficiency of CO2 trapping in clathrates and 
sis its consecutive impoverishment in Lake Vostok. 

519 The comparison of our calculated CO2 and CH4 mole 

520 fractions in Lake Vostok with future in situ measure- 

521 ments will allow to discriminate between their different 

522 supply sources. Indeed CH4 and CO2 abundances in 

523 Lake Vostok will vary according to the amount and type 

524 of biological and abiotic activity affecting the speciation 

525 of carbon. If biological activity and if water-rock inter- 

526 actions are minimal, then the concentration of both of 

527 these gases should approach the amount present in the 

528 gases at the time of Lake Vostok's formation. However, 

529 microbes, if present (as is highly likely - see Priscu et 

530 al. 1999; Karl et al. 1999), would significantly alter the 

531 abundances of these gases in the fluids. If water-rock 

532 interactions are minimal, then microbes should likely 

533 consume methane and other reducing compounds, ox- 

534 idizing them to form CO2. As such, the CO2 content 

535 of the fluid should increase relative to the amount pre- 
sse dieted, and the CH 4 content should decrease, especially 
537 because O2 would be abundant in these oxidizing wa- 
sss ters. However, if water-rock interactions are present 

539 within the lake, then chemolithotrophic microbes might 

540 survive in these environments, as suggested by Bulat et 

541 al. (2004), and may perhaps accelerate oxidation of rock 

542 material with dissolved O2. The ability of microbes to 

543 do so would be contingent on the composition of the 

544 rocks at the base of the lake. 

545 Lake Vostok may be compared favorably with Europa 

546 in terms of pressure, temperature, and potentially wa- 

547 ter composition. The pressure at the water-ice layer 

548 of Lake Vostok would be comparable to a water-ice 

549 layer on Europa occurring at 30 km. Furthermore, 

550 the presently calculated quantity of O2 dissolved within 

551 Lake Vostok is also similar to those predicted by Hand 

552 et al. (2006) and Greenberg (2010) on Europa, in which 

553 O2 would be an abundant constituent of the water. Al- 

554 though the composition of gases present on Europa is 

555 unclear, O2 has been detected on the surface of Eu- 

556 ropa (Hall et al. 1995; Carlson et al. 1999). This 2 

557 is likely transported from the surface to the subsurface 



558 ocean. The gas composition of other species is less cer- 

559 tain, however it is also highly likely that the original 

560 icy material from which Europa formed also had sev- 

561 eral gas clathrates, including those of the noble gases 

562 (Kereszturia & Keszthelyib 2012; Mousis & Gautier 

563 2004; Mousis & Alibert 2006; Mousis et al. 2011). 

564 These gases would be unlikely to have escaped from Eu- 

565 ropa, and hence should still be present in some quantity 

566 in the subsurface ocean. As the water-rock interactions 

567 at the base of the lake are unlikely to produce abundant 

568 hydrothermal systems as the underlying rock is felsic or 

569 silicic in composition (Van der Fliert et al. 2008), there 

570 is little likelihood of significant acidification occurring 

571 as O2 reacts with H2S produced by hydrothermal sys- 

572 tems, as suggested for Europa by Pasek & Greenberg 

573 (2012). For these reasons, water in Lake Vostok and in 

574 Europa is likely close to neutral in pH and highly oxi- 

575 dizing. Because of the similitude between Lake Vostok 

576 and Europa, our approach would thus be a very useful 

577 tool for accurate predictions of the composition of Eu- 

578 ropa's internal ocean, if validated by comparison with 

579 in situ measurements in Lake Vostok. 
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